Phase diagram for interacting Bose gases 
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We propose a new form of the inversion method in terms of a selfenergy expansion to access 
the phase diagram of the Bose-Einstein transition. The dependence of the critical temperature 
on the interaction parameter is calculated. This is discussed with the help of a new condition for 
Bose-Einstein condensation in interacting systems which follows from the pole of the T-matrix in 
the same way as from the divergence of the medium-dependent scattering length. A many-body 
approximation consisting of screened ladder diagrams is proposed which describes the Monte Carlo 
data more appropriately. The specific results are that a non-selfconsistent T-matrix leads to a linear 
coefficient in leading order of 4.7, the screened ladder approximation to 2.3, and the selfconsistent 
T-matrix due to the effective mass to a coefficient of 1.3 close to the Monte Carlo data. 

PACS numbers: 03.75.Hh, 05.30.Jp, 05.30.-d, 12.38.Cy,64.10.+h 



I. INTRODUCTION 

Interacting Bose gases have become a topic of current 
interest triggered by the experimental table-top demon- 
stration of Bose condensations [HQ. Especially, it is of 
great importance to know the behavior of the critical 
temperature in dependence on the interaction strength 
and the density. The fluctuations beyond the mean field 
are especially significant since they establish the devia- 
tions of the critical temperature from the noninteracting 
value and might be observable by adiabatically convert- 
ing down the trapping frequency |3J] . The measurement of 
the critical temperature of a trapped weakly interacting 
87 Rb gas differs from the ideal gas value already by two 
standard deviations Q . Therefore it is of great interest 
to understand this change even in bulk Bose systems. 

The density no of the noninteracting bosons with mass 
m and temperature T reads 



,'rnk B T 
no(T) = ( ^V 9{e) 



3/2 
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where g(x) = |[P 3 / 2 (x)] 2 ' 3 is given in terms of the poly- 



log function P%/2{x 
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and we denote the 



negative ratio of the chemical potential to the tempera- 
ture by e = —fi/ksT. The critical temperature T)(n) is 
given by ^B) for e — > 0. 

The critical temperature, T c , of the Bose system with 
repulsive interaction is deviating from the free one, To, 
and is a function of the free scattering length clq times 
the third root of the density. One usually discusses this 
deviation on condition that the density of the interacting 
system at the critical temperature, n(ao,T c ), should be 
equal to the density of the noninteracting Bose system 
no (To) at the critical temperature To. According to JT]) 
we have n (T c ) = (T c /T ) 3/2 n Q (T ) and one obtains @ 



Tc 
To 



n a (T c ) 
n{a ,T c ) 



2/3 



(2) 



The change of the critical temperature (T c — T n )/T 
in dependence on the coupling strength yo = aon 1 / 3 has 
been found to be linear in leading order. Monte-Carlo 
(MC) simulations and theoretical calculations by many 
groups have confirmed 

^-^« Cl yo + e% 2 )- (3) 
-to 

The sometimes occurring square root dependence has 
turned out to be an artifact of the first order virial expan- 
sion @. Although there is now an agreement about the 
linear slope, the actual value of c\ varies from 0.34 to 3.8 
as discussed in The situation is different for finite 

Bose systems contained in a trap. The increase of inter- 
action lowers the density in the center of the trap and 
reduces effectively the critical temperature [4J such that 
Ci = —0.93 was found in [9J. Here we want to consider 
only bulk Bose gases. 

Let us shortly sketch the different expansion schemes. 
The convergence of such expansions is discussed in 10] in 
terms of the large- N expansion of scalar field theory. The 
optimized linear expansion method [llj yielded c\ = 3.06 
for two- loop contributions and ci = 1.3 for six loops [i~2| . 
Though this (S-expansion works well in quantum mechan- 
ics it fails to converge in quantum field theory [131 ] . The 
l/N expansion [14] provides c\ = 2.33. This seems to 
be in agreement with the non-selfconsistent summation 
of bubble diagrams [y] which gives ci = ^£(3/2) ~ 4 / 3 « 
2.33. The same result has also been found with a T- 
matrix approximation [lj| . Older MC data [l6| show 
similar values such as c\ — 2.30 ± 0.25 while newer ones, 

din, 

report c\ ps 1.3. The latter result can also be 
found by variational perturbation theory [l!| yielding 
Ci = 1.23 ± 0.12. The same coefficient has been obtained 
by exact renormalization group calculations (20ll2~H which 
provide the momentum dependence of the self energy at 
zero frequency as well. The usage of Ursell operators [22j 
has given an even smaller value c\ — 0.7. 

Though the weak coupling behavior of the critical tem- 
perature can be considered as settled there is consider- 
ably less known about the strong coupling behavior. In 
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[23| a nonlinear behavior of the critical temperature was 
proposed; first the interaction leads to an increase of the 
critical temperature up to a maximal one. A further 
increase of interaction strength or correlations then de- 
creases the critical temperature until the Bose conden- 
sation vanishes at a maximal coupling parameter y§ rlt . 
One finds the interesting physics that too strong repul- 
sion prevents the Bose condensation. 



The variational perturbation method developed in [241 ] 
and advocated in (23| allows one to sum up higher or- 
der diagrams in order to account for strong coupling ef- 
fects. In this paper we present an alternative method 
well known as inversion method to obtain from the lad- 
der summation diagrams a higher order perturbation re- 
sult by inversion. The motivation to propose yet an- 
other theoretical method besides the already mentioned 
Monte-Carlo simulations and the variational perturba- 
tion method is twofold. First the strong coupling limit is 
a problem of theoretical challenge which is not yet solved. 
Different schemes favor different summations of diagrams 
and therefore weight different processes differently We 
think that different paths and comparing the results be- 
tween them might be enlightening for understanding the 
physical mechanisms. The second motivation is a pure 
methodological one. The proposed inversion method has 
led already to remarkable results in other fields. There- 
fore we think it is interesting to try how far this method 
can be applied fruitfully to the highly nontrivial problem 
of interacting Bose gases. We expect from this method 
theoretical insights into the systematics of higher order 
processes which have to be considered for the descrip- 
tion of strongly coupled systems showing spontaneous 
symmetry breaking. Here the inversion method seems 
to have its advantages since it allows explicitly for sym- 
metry breaking. This is, however, beyond the scope of 
the present work. 



Below we discuss first the method to extract the phase 
diagram within the standard ladder approximation solv- 
ing the Feynman-Galitzky form of the Bethe-Salpeter 
equation. This will be performed with the help of the 
separable potential which is then used in the limit of the 
contact interaction. This way is mathematically conve- 
nient since we can renormalize the coupling strength to 
the physical scattering length avoiding divergences. Then 
we apply the inversion method starting from the ladder 
approximation. This leads us to a twofold limit of the 
resulting series, one of which shows the nontrivial be- 
havior in the phase diagram due to symmetry breaking. 
We compare this result with the Monte Carlo data and 
other theoretical approaches in the literature. In the last 
chapter we discuss improvements, mainly the summation 
of ring diagrams connected with ladder diagrams which 
will modify the linear slope but will not change the scaled 
nonlinear curve. Possible further improvements are out- 
lined in the summary. 



II. INVERSION METHOD 

The inversion method we use has been discussed thor- 
oughly in the centennial review paper of Fukuda et al. 
[25| |. We consider a system which can undergo a symme- 
try breaking due to a small external perturbation such 
that a phase transition is possible. The corresponding 
order parameter <f> shall depend on an external source 
parameter J. We can expand the order parameter in 
terms of the coupling g to the source via 



/(■/) = £<?"/»(./) 



(4) 



n=0 



which corresponds to a perturbation series. Since we can 
perform only a calculation to finite order in g, a vanish- 
ing J means also a vanishing order parameter <p and we 
are not able to describe the phase transition as a proper 
spontaneously broken symmetry. The situation can be 
improved if (U]) is inverted 



(•5) 



For any expansion with a finite number of /„ we can 
calculate the inverted coefficients h n by introducing (Et 
into and comparing the powers of g to obtain [25| 



h (4>) 

hi(4>) 



fo(<t>) 

h(J) 



etc. 



(6) 



J=h (<P) 



If we now set J = in the inverted series, the nontriv- 
ial solution (f> ^ is obtainable though with finite-order 
calculation of the h n in ([5]). The finite truncation in 
the inverted series (JSJ corresponds to an infinite series 
summation in For further details and successful ap- 
plications of this method see (25| . 

We will use the ladder summation to establish a sum 
for the scattering length and will invert this series with 
respect to the selfenergy in order to obtain a nonlinear 
phase transition curve, i.e. the dependence of the critical 
temperature on the coupling strength. 



III. MANY-BODY T-MATRIX 
APPROXIMATION AND BETHE-SALPETER 
EQUATION 

The T-matrix depends on the frequency w/h and the 
center-of-mass momentum K and reads according to fig- 
ure [1] 

TV (w, K) = V vv , +Y, v pi ]_ +f ^ f ^^ T iP' ( w > K ) 



(7) 



where we indicate the difference momenta of the incom- 
ing and the outgoing channel by subscript, the energy 
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FIG. 1: The diagrams of the T-matrix approximation where 
the difference momenta are indicated as well as the momenta 
of the two intermediate propagators. 



dispersion is assumed to be e p — p 2 /2m and the Bose 
functions are denoted by fk- 

This Bethe-Salpeter equation can be solved by a sepa- 
rable potential V pp i = \g p g p > with a coupling constant A 
and a Yamaguchi form factor g p = l/(p 2 + (3 2 ) where /3 
describes the finite range of the potential. One obtains 



The scattering phase 
ImT 



Xg P g P < 



q 



\J{uj,K) 
! + /«- 



+ /*- 



Ejc. 



£k 



+q 



IT] 



tan S = 



ReT 



p = p\K = 



ak + 0{k 2 ) 



(8) 



(9) 



is linked to the scattering length by the small momentum 
expansion, p = hk. From ([8]) we obtain 26] 



+ 2/o) 



A jf dpjjp^yj 



(10) 



Now we can perform the limit to the contact potential 
using a /3-dependent coupling constant A(/3) in such a 
way that the free two-particle scattering length do in the 
absence of many-body effects, fk — * 0, is reproduced [26| . 
From ([Toll we have 



2ir 2 h 3 



A(a ,/3) 



/ dp { p'il p -i)i 



(11) 



which we use to eliminate A in (fTTj)) to obtain o(oo,/3). 
Performing now the limit (3 — > oo we obtain a well de- 
fined expression for the scattering length with many body 
effects 



a 
Oo 



l + 2/o 



o 



(12) 



This scattering length depends on the chemical potential 
and on the temperature, a(/i, T) via the Bose function f p . 
This renormalization procedure is completely equivalent 



to a momentum cut-off used in [27]. Higher order par- 
tial wave contact interactions have been discussed in [28[ 
which are important for constructing pseudopotentials 
[29l ] . For further relations of the T-matrices in different 
dimensions see (30| . 

In [26| the expression (fl~2)) was discussed and it was 
found that this medium-dependent scattering length di- 
verges at the Bose-Einstein transition characterized by 
the critical density and temperature. We will employ 
(fl"2)) as a good measure of approaching the Bose-Einstein 
condensate. With the help of ([T]) we can rewrite (fT2")l into 
an equation for the dimensionless interaction parameter 
= an 1 / 3 versus the free one yo = aon^ 3 , 



JJ_ 

yo 



coth(^) 
l + 4y ff'(es) 



(13) 



on the condition n(ao,T c ) — no(To). The ratio of the 
temperatures x = T c /Tq for interacting to noninteracting 
Bose systems follows from (JT]) and (|2j as 



9(0) 



(14) 



where = e + S(es, yo)/^B^ denotes the functional 
dependence of the quasiparticle energy in the argument 
of the distribution function on the selfenergy. 



IV. BOSE-EINSTEIN CONDENSATION 
BORDER LINE 

A. T-matrix approximation 

Our aim is to extract the phase diagram of do in terms 
of x = T c /Tq. A simple elimination of e in the system 
of equations (fl"3"|) and (fTI)) would not lead to an answer 
since the ladder approximation is a weak coupling theory 
[3ll ] and cannot lead to the proper description of strong 
coupling. The simple elimination e = e(y) — > x = x[e(jj)] 
would correspond to the direct series ((4]). We will present 
how the strong coupling limit can be reached in the 
next chapter. Here we derive a proper condition for the 
Bose-Einstein condensation border line. This will be per- 
formed in two ways, first by the observation of a diverg- 
ing medium dependent scattering length and second by 
the condition of the appearance of a Bose pole in the T- 
matrix. Both conditions will lead to the same criterion. 

Starting we take into account the condition of being 
as near as possible to the phase separation line. We use 
the fact that the medium-dependent scattering length di- 
verges at this phase separation line such that we obtain 
for critical eff ' -> e from (JTSj) and (M]) 



yo 



y 



1 



coth^ - 4y ff '(e E ) 



Vo 



5(0) 
9(e) 



(15) 
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as the system of equations which gives the separa- 
tion line in the phase diagram. The meaning of e = 
—/i(ao,T c )/kBT c + Yj/ksTc is the critical chemical po- 
tential and the selfenergy for the interacting Bose gas 
divided by the critical temperature. 

The same condition (fT5|) for the formation of a conden- 
sate can be found by calculating the binding pole of the 
T- matrix ([5]) . With the renormalization (fTTj) the bosonic 
T-matrix (151 for contact interaction reads 



T pp ,(w,Q) = - 



4irh 2 a 



1 



TO l + J(w,Q) 
In the center of mass frame one finds 



(16) 



OO 



q 2 f q 
q 2 — muj 



clq J y 7 — Lum . uj < 



(17) 



In the dilute limit, f q — > 0, the free two-particle phase 
shift of © starts from zero for positive ao and from tt 
for negative scattering length indicating a bound state 
according to the Levinson theorem [32|. For attractive 
interactions, ao < 0, we have therefore a bound state 
pole at J(—Eb,0) = —1 with the binding energy Eb = 
H 2 /maQ. With the sign of the scattering length in Q 
we follow the convention of [13] which has an opposite 
sign compared to [33| . Generally odd numbers of bound 
states correspond to negative scattering lengths and even 
numbers to positive scattering lengths. 

We are interested here in the Bose condensation, that 
means in a pole for repulsive interaction, ao > 0. In this 
case there is no bound state, but a Bose-condensation 
pole possible at J(—Eb,0) = —1. We can estimate 

OO 

J( W <0,0) > J(0,0) = -^ j dqf q =Ay a g'{e). (18) 



The T-matrix (fTrJj) has therefore a Bose-condensation 
pole if J(0, 0) < — 1 which means in view of l|T5|) that 



2/o > 



1 



(19) 



This is the same condition as (TT51) and one can see that 
the existence of a pole in the T-matrix is equivalent to 
the existence of a Bose condensation. In other words 
this pole represents the medium dependence of the Bose 
condensation in T-matrix approximation. 

Let us shortly discuss the physical content of the pole 
condition of the T-matrix. The result according to (JT7J) 
is plotted in figure [5] versus the dimensionless interaction 
parameter. For e = — /i/fc^T — * the curve starts at zero 
which means that any small interaction leads to a pos- 
itive Bose condensation energy while for noninteracting 
systems the condensation energy remains at zero. For the 




FIG. 2: The condition for appearance of a Bose-condensation 
pole, J(-E b ,0) = -1 according to l[T7|). The Bose- 
condensation energy is plotted versus the dimensionless in- 
teraction parameter for different ratios of the chemical poten- 
tial to the temperature, e = — /n/fesT. The inset shows the 
minimal interaction parameter as a function of e. 



critical e = the curve approaches Eb = 1.011 x 2ksT c 
for large interactions which shows that a further increase 
of the Bose condensation energy is not possible. This 
gives already a hint that a too strong interaction is not 
increasing Bose condensation further. In fact we will find 
in the next chapter within the strong coupling summation 
that the Bose condensation will decrease with increasing 
coupling parameter. Here in the weak coupling limit it 
stays constant. 

A further interesting observation is that for lower den- 
sities e > there appears a minimal coupling parameter 
above which only a positive Bose condensation energy is 
possible. Below this minimal interaction parameter we 
have a strong imaginary part of the pole (|17p indicat- 
ing an instable state with a finite lifetime. This minimal 
interaction parameter as a function of e is seen in the 
inset of figure O In fact it rises exponentially such that 
for a given interaction strength we have an upper limit 
in e = —fJ,/T up to which Bose-Einstein condensation is 
possible. 



B. Inverse selfenergy expansion 

We can describe the strong coupling limit by the proper 
inversion method as demonstrated in the following. Anal- 
ogously to the inversion method discussed in chapter II 
we identify now the order parameter (j) with x, the pa- 
rameter J with e and the expansion parameter g with 
Vo- 

The separation line of Bose-Einstein condensation ap- 
pears for vanishing e. Inverting the small e expansion of 
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variational perturbation method [H( (thin dashed line in 
figure [3]) is nearer to the MC data. This will be improved 
in the next chapter. 

Let us stress again that this inversion procedure is not 
a simple expansion of x(y§) in powers of ?/§ but primar- 
ily an expansion in e(x). In fact we had to expand up to 
second order in e in order to obtain the third-order coeffi- 
cient in (|20|) , to third order in e in order to obtain fourth 
order in (|20p. etc. This is the reason why we do not 
get a logarithmic term for the second term as obtained 
perturbatively in second order in (3?| . 



FURTHER IMPROVEMENTS 



A. Towards selfconsistency 



FIG. 3: The phase diagram of the scaled dimensionless scat- 
tering length versus the critical temperature. The MC data 
are from 34,351 and the variational improved perturbation 
theory from [8ll23l|. 



the function x(e) in the second line of (fl"5f the resulting 
e(x) expansion corresponds to the finite series In- 
serting this into the first line of Q15p gives ?/o( x ) which 
again is inverted leading to 



x = 1 + c lV c + c 2 (y c Q y + c 3 {y c Q y + . 



(20) 



with 

ci = 

C2 = 
C3 = 



16tt 



3C(1) 4/3 



4.657 



8.325 



327T (14tt + 9C(|)C(|)) 
9C(f) 8/3 

512^(88^ + 108K(|)C(|)+27C(i) 2 C(|) 2 



81C( 



= -14.032 



(21) 



where we have derived the expansion up to 8th order 
in e in order to achieve sufficient convergence. It has 
to be remarked that the series possesses diverging terms 
with alternating signs which nevertheless converge to two 
limiting curves. 

The results (|20|) are plotted in figure [3] One of the two 
limiting curves (dotted line), called 2 nd solution, corre- 
sponds to the naive elimination x = x[e(?;g)] discussed 
above. The thick dashed line shows the new result due to 
the inversion method which goes beyond the level of the 
T-matrix summation of diagrams. We see that the tran- 
sition curve bends back and a maximal scattering length 
occurs below which Bose condensation can happen. It 
is instructive to compare also the absolute values of the 
maximal critical parameters where the Bose condensa- 
tion vanishes. These are presented in table |TJ Our maxi- 
mal critical temperature is too high. In this respect the 



Though the figure [3J describes the data quite well it 
is mainly due to the scaling of the a axis to one. It 
is more instructive to compare the leading order ci ex- 
plicitly. The non-selfconsistent T-matrix leads to an in- 
crease which is too large by nearly a factor 2 compared 

with @0iiii6i. We can then improve this result by the 

selfconsistent T-matrix. This is achieved if we replace 
the free dispersion p 2 /2m in the arguments of the Bose 
function in (P),© and (TTJ)), by the quasi-particle one e p 
which is a solution of the selfconsistent equation 



2m 



(22) 



The retarded selfenergy reads [37 



£(w,fc) 



Tk- P k-p {lo + e p , k + p)f p 



E 



\Tk-p k^p_ (e k -q+£p+q,k+p)\ 2 fk-qfp+q 



£k-q + £p+q - e p -LJ + ir) 



(23) 



with the retarded renormalized T-matrix ([5]) in the limit 
of contact interaction 



Tk-p k-p (uj, k+p) = 

2 ' 2 y 



4irh 2 aQ/m 



l+47rft 2 aoE„< 4t 2 

u \ p' 2 p' 2 -\-EZ- 



(24) 



Since the medium-dependent scattering length (|12p is de- 
termined by the small k expansion and the poles of the 
different occurring Bose functions appear for small e and 
consequently small momenta, we consider the quasipar- 
ticle features only at small momenta. One should note 
that we are still in the gas-like limit for the quasiparti- 
cle dispersion as can be seen from its quadratic behavior. 
Further improvements require the consideration of the 
condensate explicitly [38| leading to linear phonon-like 
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(Tc/Ti) )max 


cut 

yo 


Cl 


T-matrix 


1.56 


0.16 


4.66 


T- matrix (m*(e)) 
T-matrix (m*(0)) 


1.22 


0.16 


1.32 


1.07 


0.16 


-2.13 


Screened T-matrix 


1.56 


0.32 


2.33 


MC data [34] 


1.06 




0.34±0.08 


MC data [16] 






2.30±0.25 


MC data [17] 






1.32±0.02 


MC data [18] 






1.29±0.05 



TABLE I: The maximal critical temperature, the critical cou- 
pling parameter where the Bose condensation vanishes and 
the linear coefficient from the critical temperature according 
to figure [3] for the three different approximations used in the 
text. 



For our present discussion we have used only the first 
part of the selfenergy (f2"3")l . The system of equations ([55]) 
is solved and the resulting curve approaches the MC data 
as can be seen in figure[3] (thick solid line). Especially the 
maximal critical temperature is nearer to the MC data as 
can be seen in tableQ] Now the effective mass depends on 
e. This leads to a delicate numerical balance as illustrated 
by the curve for the effective mass at fixed e = which 
gives the best overall agreement with the data points 
but shows a negative c\ . Calculations with higher-order 
approximation of the selfenergy are therefore certainly 
necessary to achieve a sufficient convergence. 

B. Screened ladder diagrams 



disp ersions [39j up to roton minima at higher momenta 
[4fJ ] . For the present level we restrict our investigation to 
the quasiparticle renormalization at small momenta and 
quadratic dispersion. 

The one-particle dispersion relation takes the form 



2m 



Re£(0,0) 



p 2 d 

u 2m 



ReS(w',p') 



o( P a 



p'—u'—Q 



2m* 



aReS(0,0) + O(p 3 



with 



1 - #7Re£(u/,0)| , . 



(25) 



(26) 



and the effective mass 



d 



m*(e,y ) 



a 1+ Re£(a/,p') 



QE1. 



d d 



\ 2m* 



7 ReE(w',j/) 



to'—p'—O 



(27) 



The second expression is the selfconsistent one taking 
into account that the effective mass appears in all ex- 
pressions of the selfenergy. The constant selfenergy shift 
aReE(0,0) can be included into the definition of the 
chemical potential in the same way as the mean-field con- 
tribution which does not lead to any change in the critical 
temperature. All what is left from selfconsistency is the 
effective mass. 

It is not difficult to see from ([T]) and @ that in the 
second line of (fl~5|) an additional factor m/m* has to ap- 
pear 



.i/o 



W) 

g(0) m 
g(e) m*(e,yg) 



(28) 



We can now further improve the many-body approxi- 
mation by using a ring summation of diagrams for the po- 
tential of the T-matrix. This corresponds to the screened 
ladder approximation such that the bare potential in the 
T-matrix is replaced by a screened one [4l|. In order to 
maintain separability of the T-matrix we consider this 
screened potential in separable form as a result of the 
screening 



V pp > — V pp > V f 



/Pi— f P - 

J 2+p 2 -p y (nq\ 

pp . v pp' \^- J ) 
ep+p- ££-p- w + ir] 



which is solved yielding 



y- _ ^9p9p> 

pp ' i + xn(u, P ) 



n(o;,p) - J2 



en 



+p~ EE^p - UJ + 17] 



(30) 



We see that the screening corresponds to a scaling of the 
coupling constant in the T-matrix ((7|) 



A 



A 



l + XU(uj,p) 



(31) 



and the same solution as ([S]) appears but with a new 
momentum- and energy-dependent coupling. We can 
now follow the same procedure as outlined above in using 
a /3-dependent coupling constant to perform the contact 
potential limit (3 — > oo. For this aim we need the small 
momentum expansion 



dpf P +^+o(p 2 ). 



87Tft 3 /? 4 



The resulting formula analogous to (fT2"]l reads now 

l + 3/o 



a 

do 



o 



(32) 



(33) 
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which leads to 



Vo = 



2<?'(e) 
.9(0) 



9% 



(34) 



instead of (fl~5|) . Therefore an additional factor 1/2 ap- 
pears in the linear coefficient (f2"Tj) . 



8- 



3(( 



,4/3 



= 2.33, 



(35) 



a result in remarkably good agreement with [6 HT^fl5| and 
the MC data [16|. The critical values are found in ta- 
ble HI The screened ladder approximation leads to twice 
the maximal critical coupling strength but the same max- 
imal critical temperature compared to the T-matrix ap- 
proximation. We note that further improvements should 
be possible again by considering the selfconsistency, i.e. 
the effective mass. 



VI. SUMMARY 

We have discussed interacting Bose systems with a re- 
pulsive interaction. The appearance of the Bose-Einstein 



condensation becomes dependent on the interaction pa- 
rameter. We derive the same condition for Bose-Einstein 
condensation both from the pole of the T-matrix and 
from the observation that the medium-dependent scat- 
tering length is diverging at the Bose condensation. 
This condition is used to produce a perturbation series 
of higher order than the original T-matrix summation 
adopting the well known method of inversion. With the 
help of this method we obtain a strong coupling result out 
of the weak coupling T-matrix summation. We are able 
to describe the nonlinear behavior of the critical temper- 
ature. The critical temperature rises linearly with the 
interaction parameter in leading order reaching a maxi- 
mal value and is reduced for higher interaction strength. 
At an upper specific interaction parameter Bose conden- 
sation is not possible anymore. The system is too repul- 
sively correlated for allowing condensation. 

We improve the original weak-coupling series further 
by using screened vertexes which reproduce the leading 
order coefficient of the MC data. Further improvements 
towards the nonlinear curve of the MC data and the vari- 
ational perturbation theory are possible by the selfcon- 
sistency. Here we present as a first step an effective mass 
calculation inside the T-matrix and show that the in- 
version method leads then to a better description of the 
data. 
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